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ABSTRACT 

Compared to starburst galaxies, normal star forming galaxies have been shown to display a much 
larger dispersion of the dust attenuation at fixed reddening through studies of the IRX-/3 diagram 
(the IR/UV ratio "IRX" versus the UV color "/?"). To investigate the causes of this larger dispersion 
and attempt to isolate second parameters, we have used GALEX UV, ground-based optical, and 
Spitzer infrared imaging of 8 nearby galaxies, and examined the properties of individual UV and 
24 fiui selected star forming regions. We concentrated on star-forming regions, in order to isolate 
simpler star formation histories than those that characterize whole galaxies. We find that 1) the 
dispersion is not correlated with the mean age of the stellar populations, 2) a range of dust geometries 
and dust extinction curves are the most likely causes for the observed dispersion in the IRX-beta 
diagram 3) together with some potential dilution of the most recent star-forming population by older 
unrelated bursts, at least in the case of star-forming regions within galaxies, 4) we also recover some 
general characteristics of the regions, including a tight positive correlation between the amount of 
dust attenuation and the metal content. Although generalizing our results to whole galaxies may not 
be immediate, the possibility of a range of dust extinction laws and geometries should be accounted 
for in the latter systems as well. 

Subject headings: 



1. INTRODUCTION 

Star formation is a fundamental process that governs 
the evolution of the baryonic matter in galaxies. It 
can deeply affect the host galaxy, its properties and ap- 
pearance. The newly formed populations change the 
galaxy's spectral energy distribution (SED); the met- 
als formed in the most massive stars pollute the in- 
terstellar medium and, in the presence of mechanical 
feedback from episodes of intense star formation, can 
be ejected from the galaxy and pollute the intergalac- 
tic medium. It is important to accurately measure star 
formation both as a function of spatial position within 
a galaxy and as a function of a galaxy's temporal evolu- 
tion. The cosmi c star formation ra t e density as a func- 
tion of redshift (|Madau et alj 119961 : iHopkins fc Beacoml 
120061 ) has become a classic comparison benchmark, and 
a challenge, for both semi-analytic and hydrodynami- 
cal models of structure formation and evolution (e.g. 



Katz et al.ll996UHopkins et alJ2005lStinson et~aT1l2006l : 
Cattaneo et al.ll2007l ). The youngest stellar populations, 
whose bolometric energy output is dominated by short- 
lived massive stars, emit the bulk of their energy in the 
restframe ultraviolet (UV) . Since the pioneer work of 
iDonas fc Deharvengl (|1984f ) the UV has become one of 
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the main star for mation rate (SFR ) tracers used in mod- 
ern Astronomy (jKennicuttl 119981 ). This regime, which 
has been previously difficult to observe in the local Uni- 
verse, becomes the wavelength region "par excellence" 
to investigate star formation in galaxies at intermedi- 
ate and high redshift. At these distances, the restframe 
UV is shifted to optical and near- infrared wavelengths, 
where current detectors afford the highest performance 
in terms of both sensitivity and angular resolution. Thus, 
it is no surprise that so far galaxy investigations in the 
redshift range z^ 2-7 have been dominated by restframe 
UV observations jG iavalisc o et alj 120041; iBouwens et al.1 



l2TMlMobasher et al.H2005l : ISawicki fc Thompsonll2006h . 
UV observations are the main canvass from which our 
current understanding of cosmic star formation is built. 

For over a decade now, studies of galaxies at cosmo- 
logical distanc es have employed the "st a rburst atten- 
uatio n curve" (jCalzetti et alj 11994 l2000t iMeurer et all 
119991 ) - see next paragraph for a description of the 
physical con ditions leading to su ch a law and equa- 
tion 4 from iCalzetti et al.l (|2000f ) - to correct their 
restframe UV measurements for the effects of dust ex- 
tinct i on (e.g. iSteidel et all 119961 Il999t iGiavalisco et all 
[200l iDaddi et all 120051 . 120071 . to name a few). This 
curve is a powerful tool, in that it only requires infor- 
mation on the UV colors (or UV spectral slope) of a 
galaxy to provide an extinction correction, and recover 
the intrinsic UV flux. Assuming this curve, the typical 
UV extinction correction at reds hift 3 is a factor ~ 5 
in luminosity ([Steidel et al.|[l999T) . However, in recent 
years, evidence has been accumulating that the star- 
burs t attenuation curve is not applicable to all galax- 
ies (|Buat et al.1 [2001 [20051 : iNoll et all l200l and that 
the UV extincti on is dependent on the age of the stel- 
lar populations (jCortese et alj [20081 ). More quiescently 
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star-forming galaxies show characteristics in their UV 
SED that do not lend themselves to the same extinc- 
tion "treatment" as more active starbursts. Applying 
the starburst attenuation curve to these more quiescent 
galaxies leads to up to an order of magnitude overesti- 
mate of the UV luminosity (and of the SFR) . This result 
clearly represents a problem as deeper observations probe 
past the most powerful starburst to more normal galaxies 
with relatively lower SFR at increasingly higher redshifts. 
Along these lines, current determinations of UV lumi- 
nosity functions at cosmological distances may be viewed 
with some suspicion, as all of them need some form of ex- 
tinct i on correction (iSteidel et al.l 119991; iG iavalisco et al.l 
[200l [Sawicki fc Thompsonll2006l ; iTeplitz et al.ll2OO60 . 

One of the prescriptions of the starburst attenuation 
curve is based on a very simple, fully-empirical result. 
The ratio L(IR)/L(UV) is strongly correlated with the 
measured UV spectral slope (3 (or the UV color) in lo- 
cal, UV-bright star burst galaxies (|Meurer et al.l 119991 ; 
iCalzetti et atll2000t) . The L(IR)/L(UV) ratio is a mea- 
sure of the total dust opacity that affects the UV emis- 
sion in a galaxy, or more generally, in a stellar system; 
the UV light absorbed by dust is re-emitted in the IR, 
independently of the details of the star-dust distribu- 
tion, the dust properties, the mctallicity, the stellar ini- 
tial mas s function, or the star formation hist ory of the 
system (| Gordon et al.ll2000t iBuat et all 120051 ). Meurer 
et al.'s result, commonly called the IRX-/? relation, to- 
gether with a host of previous results that showed strong 
correlations between (3 and a number of dust redden- 
ing measures, form the basis for an extinction correc- 
tion prescription that obviously a pplies to local (and 
possibly distant) starburst galaxies (jCalzetti et al.|[l99"4l 
119961: iCalzettil fl997l 120011 ). The properties of the star- 
burst attenuation curve imply specific conditions for the 
distri bution of dust and stars within the starburst re- 
gion (ICalzetti et al l [l99i [19961 iWitt fc Gordon! 120001: 
iCharlot fc Fallll200d ICalzettill200lD . The geometry that 
can account for all observed properties at UV, optical, 
near-IR, and far-IR wavelengths is that of a dust distribu- 
tion foreground to the stellar population responsible for 
the bulk of the UV emission. The term foreground is here 
loosely used to indicate any geometry that puts the stars 
in the center of a shell-like dust distribution, and does 
not make any inference on the dust distribution (which 
can be homogeneous or clumpy). Starbursts provide a 
mechanism for creating such specific dust-star geometry: 
massive star wind and supernova feedback can create ex- 
panding bubbles of hot gas that sweep away a significant 
fraction of the dust and gas that resides in the starburst 
site itself. Within this scenario, the gas and dust end 
up at the "edges" of the starburst site thus creating an 
inhomogencously distributed "shell" . In addition, the 
starbust extinction curve also requires a specific type of 
dust grains lacking the 2 175 A b ump dGordon et"aTI|1997f) 
like that found in th e SMC bar (jGordon fc Clavtonlll998l : 
iGordon et all[2003l ). 

The starburst attenuation curve may thus not apply 
to a generic galaxy, as the coherent feedback mechanism 
present in starbursts may be missing in other galax- 
ies. A number of studies have confirmed that more 
quiescent star-forming galaxies, although they gener- 
ally trace a correlation in the IRX-/3 plane, are also 



up to an order of magnitude less obscured than star- 
bursts for similar (3 values, and can have a larger dis- 
persion, with a facto r up to ~5 larger spread around 
the mean correlation (IBuat et all 12001 120051: iBelll 1 2002 : 
iKong et al.112001: IGordon et al.ll2 004; Sci bert et al| 2005: 
Calzetti et al.ll2005l: iBoissier et al.ll2007HDale et al.ll2007 : 
Munoz-Mateos et al.l 1200*91 )! UltraLuminous Infrared 



Galaxies (ULI RGs) also deviate from the starburst IRX- 
(3 correlation (jGoldader et alj [2002), and are located 
above it, i.e., opposite to the star-forming galaxies. 
ULIRGs are dustier, more compact, and denser sys- 
tems than the UV-bright starbursts; therefore, their loca- 
tion in the IRX-/3 diagram is plausibly account ed for by 
an "h omogeneous" mixture of dust and stars (|Calzettil 
120011 ). ULIRGs likely represent less of a concern for 
the use of the UV emission as a SFR tracer, as in gen- 
eral they will not be detectable at rcstframc UV wave- 
lengths in cosmological samples. Quiescent star-forming 
galaxies are a far larger concern, as they are relatively 
bright in the UV and become increasingly detected as 
surveys go deeper and, by the nature of the luminos- 
ity function, are more numerous than starbursts. Alas, 
for these galaxies, the next step in the investigation, un- 
derstanding and quantifying the nature of the deviation 
from the star bursts ' IRX- /3 correlation, has so far eluded 
any attempt. iBelll l)2002l ) suggested that a combination 
of variations in the dust geometry and the characteris- 
tics (age, etc.) of the stellar populations could be the 
reason for the deviation observed in the HII regions of 
the Large Magellanic Cloud. The sparse data available, 
however, prevented further insights into the problem. 
Simila r limitations were encountered by IGordon et al.l 
(|2004l ) in their analysis of the star-forming regions in 
M 81. IKong et al.l (f200l . using a non-homogeneous UV 
dataset, suggested that the deviations of a large sam- 
ple of quiescent star-forming galaxies from the starburst 
relation could be parametrized by the b-parameter, the 
ratio of the curre nt to past ave r age st ar formation rate. 
However, neither ISeibert et al.l (|2005l ). using the more 
homogeneous GALEX All-sk y Imaging Survey sa mple 
plus IR data from IRAS, nor I Johnson et aD 1)20071 ) . us- 
ing SDSS data, Spitzer, and spectroscopic observations, 
confirmed Kong et al.'s suggestion. Some of these works 
have also produced "mean" obscuration curves for quies- 
cent star-forming galaxies, but the intrinsic spread of the 
data around these mean curves is more than an order of 
magnitude, versus ~2-3 typical of the starbursts. This 
suggests the existenc e of a "second pa rameter" for the 
star- forming galaxies l|Kong et al.ll2004f ). which needs to 
be identified and understood if the UV is going to be 
used as a tracer of SFR. 

In this paper, we investigate the possible reason for 
the deviation from the starburst IRX-/3 relation concen- 
trating on sub-kpc star forming regions within galaxies, 
which provide simpler samples of stellar population than 
whole galaxies. We use GALEX observations to mea- 
sure the UV luminosity, Spitzer ones to derive the total 
infrared luminosity and optical ones to estimate the age. 

In section 2 we present the sample of galaxies selected 
for this study, in section 3 we present the UV, IR and 
optical observations, in section 4 we present the data pro- 
cessing and photometry methods carried out, in section 
5 we present how we model the star forming regions, in 
section 6 we proceed on to test the validity of the age 
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as the "second parameter" , in section 7 we test different 
extinction laws, in section 8 we test the metallicity, we 
discuss the results presented in section 9 and we finally 
conclude in section 10. 

2. SAMPLE OF GALAXIES 

Since we wish to investigate individual star forming 
regions within galaxies, we select galaxies where such re- 
gions can be resolved. Indeed the typical resolution in 
mid-infrared at 24 [im with Spitzer MIPS or ultraviolet 
with GALEX - which trace star formation - is about 5- 
6". The typical size of a giant HII complex, about 300 pc 
can be resolved at a distance up to 10 Mpc at this res- 
olution. Within this size at most a few stellar clusters 
may reside and contribute to the ultraviolet and infrared 
emission. Thus this constraint on the distance ensures 
that there are as few clusters as possible in an individual 
giant HII region. On the other hand we also exclude very 
nearby galaxies that would provide a very high level of 
detail. We have applied the following criteria to select 
the sample of galaxies for this study: 1) a distance closer 
than 10 Mpc, 2) to be star-forming spiral galaxies, 3) to 
cover a wide range in total SFR and SFR surface density, 

4) not to be too inclined (minor axis/major axis > 0.5), 

5) have GALEX, Spitzer data available. 
The final sample consists of 8 star- forming galaxies, 

details are given in Table [1] Their star formation rates 
range from 0.004 M Q yr _1 kpc~ 2 to 0.046 M Q yr _1 kpc" 2 
with an average of 0.020 ± 0.021 M Q yr _1 kpc -2 calcu- 
lated over the area defined by i?25- 

3. OBSERVATIONS 
3.1. GALEX 

We have obtained far and near ultraviolet processed 
count and intensity maps from the GALEX archives. 
Thanks to its large field of view of 1.24°, it can en- 
compass the full extent of the selected galaxies with a 
single pointing. The far ultraviolet band is centered on 
A e ff=151 nm whereas the near ultraviolet one is centered 
on A e -g?=227 nm. Detailed information about the instru- 
ments are provided in lMartin et al.l (|2005[ ). 

Galaxies have been observed as part of the Nearby 
Galaxy Survey (NGC 0628, NGC 2403, NGC 5055, 
NGC 5194, NGC 5457 and NGC 7793) and Guest in- 
vestigators programs (program 71 cycle 1 for NGC 3031 
and program 5 in cycle 3 for NGC 5236). 

3.2. Ground based optical 

Ground based optical images consist of narrow-band 
Ha observations as well as broad band U, B and R ob- 
servations. 

Most Ha images have been obtained from the 
Spitzer Infrared Ne arby Galaxies Survey (SINGS, 
iKennicutt et alJ l2003h archives, fifth delivery, final en- 
hanced products released on 2007-04-10. Observations 
have been carried out at the National Optical Astronomy 
Observatory, in 2001 on the Kitt Peak National Obser- 
vatory (KPNO) 2.1 m and Cerro Tololo Inter- American 
Observatory (CTIO) 1.5 m telescopes. More details arc 
available in the SINGS "User's Guide" 8 . For NGC 5236, 

8 http: //data. spitzer . caltech . edu/popular/sings/\20070410. 



the data were obtained on the Bok 2.3 m telescope on 
2005-05-10. 

Continuum emission has been subtracted from the Ha 
image using the corresponding scaled R band image fol- 
lowing the method described in Hong et al. (in prepara- 
tion). 

SINGS spectroscopic observations (Moustakas et al., in 
preparation) have been used to infer the mean [Nil] /[Ha] 
ratio within the galaxies. The measured fluxes have 
been corrected for contamination by the [Nil] lines 
accordingly. As no spectrum was available for the 
high-mctallicity galaxy NGC 5236, w e have assumed 
[NII1/ [Ha]=0.5, close to 0.53 found by IKennicutt et all 
(|2008f ) which has been estimated based on a scaling re- 
lationship between M(B) and [Nil] as described in ap- 
pendix B of that paper. Although the relative strength 
of the [Nil] lines can change from one HII region to an- 
other and can evolve as a function of the radial distance, 
we have assumed a constant ratio for each galaxy. In- 
deed very super-metallic and sub-metallic star forming 
regions - in which the relative strength of the [Nil] lines 
may deviate significantly from the mean observed values 
- are less likely to be detected in the UV for the former 
and in the mid-infrared for the latter. 

Images in the U, B and R bands have been obtained 
either from the SINGS archives - as mentioned in the pre- 
vious section - or observed in the context of the Steward 
Nearby U-band GaLa xy (SNUGL) survey, using 90prime 
(I Williams et al.ll200l . 

3.3. Spitzer 

Spitzer is a key observatory to estimate the total in- 
frared luminosity. Indeed its two imagers cover a large 
range of infrared emission, from 3.6 /im to 8.0 /jm for 
IRAC and 24 fim to 160 fxm for MIPS. The relations 
to estimate the bolometric infrared luminosity routinely 
use the IRA C 8.0 ^m and the M I PS 24 nm bands (see 
for in stance [Calzetti et al. 2005j; iPerez-Gonzalez et alJ 
120061 Boquien et al. in preparation, and references 
therein). All galaxies but two have been observed as 
part of SINGS. We use the processed legacy data pro- 
vided by the SINGS team. The only galaxies that are 
no t observed by SINGS a re NGC 5236 (also observed 
bvlEngelbracht et al.l l2008h and NGC 5457 (observed by 
IGordon et al.ll2008f ) are part of the Spitzer Local Volume 
Legacy Survey (Dale et al., 2009). 

3.4. Summary 

Images of the galaxies in FUV, Ha and 24 fj,m are 
presented in Figures [1] and O For NGC 5236 the missing 
Ha image has been substituted with the R band one. 

A summary of the bands, telescopes used, observing 
dates, exposure times and sensitivities of all observations 
are presented in Table El 

4. DATA PROCESSING AND PHOTOMETRY 

4.1. Definition of the apertures 

Given the multi-wavelength nature of this study, we 
have selected the star forming regions using the GALEX 
FUV and the Spitzer MIPS 24 /im images. Indeed, not 
only do those wavelengths trace star formation but they 
also have a similar resolution (^5-6") which ensures that 
the totality of the region is encompassed in other wave- 
_enhak»gi^iytfa^p:t)§/kiffi^ho^ as the IRAC 
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TABLE 1 
Galaxy parameters 



Name 


Morphology 




«5 a 


Distance 


PA b 


Major a 




Minor axis a 


P25 C,d 


[0/HJ d 


A[0/H] d 


[NII]/Ho 






(J2000.0) 


(J2000.0) 


(Mpc) 


(deg) 


(') 




(') 


(') 


(P = 0) 


(dox/p 25 ) 




NGC 0628 (M 74) 


SAbc 


01 36 41.8 


+ 15 47 01 


7.3 





10 


5 


9.5 


5.24 


9.16 


-0.58 


0.47 


NGC 2403 


SABcd 


07 36 51.4 


4-65 36 09 


3.13 


123 


21 


9 


12.3 


10.94 


8.88 


-0.27 


0.63 


NGC 3031 (M 81) 


SAab 


09 55 33.2 


4-69 03 55 


3.55 


155 


20 


9 


14.1 


13.46 


9.08 


-0-43 


0.28 


NGC 5055 (M 63) 


SAbc 


13 15 49.3 


4-42 01 45 


7.8 


107 


12 


6 


7.2 


6.30 


9.27 


-0.54 




NGC 5194 (M 51) 


SABbc 


13 29 52.7 


+47 11 43 


7.62 


22 


11 


2 


6.9 


5.61 


9.30 


-0.49 


0.50 


NGC 5236 (M 83) 


SABc 


13 37 00.9 


-29 51 56 


4.5 





12 


9 


11.5 


5.9 


9.24 


-0.19 


0.50 


NGC 5457 (M 101) 


SABcd 


14 03 12.6 


+54 20 57 


7.2 


39 


28 


8 


26.9 


14.1 


8.89 


-0.93 




NGC 7793 


SAd 


23 57 49.8 


-32 35 28 


3.91 


102 


9 


3 


6.3 


4.67 


8.93 


— 0.32 


0.28 



' Data from the NASA/IPAC Extragalactic Database. 
' The position angle is calculated manually adjusting an ellipst 
Radius of the major axis at the fj.^— 25 mag arcsec ^ isopho 
1 Data obtained from [Zarit sky ct "aTl |1994|) and Moustakas et j 



preparation). 
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Fig. 1. — Combined FUV (blue channel), Ha (green channel) 
and 24 (im (red channel). The top-left image in NGC 0628, the 
top-right one NGC 2403, the bottom-left one NGC 3031 and the 
bottom-right one NGC 5055. Note that the Ha image of NGC 3031 
does not cover all the galaxy and is the combination of several 
pointings which have slightly different background levels can can 
be clearly seen in the image. 

8.0 /im and the optical bands. Given that they offer a 
complementary view of star formation - far ultraviolet 
emission absorbed by the dust is recmittcd in the in- 
frared -, the selection of the star forming regions has 
been performed independently in those bands. This al- 
lows us not to miss regions that have a very low or a very 
high extinction, the former being mainly visible in FUV 
while the latter are primarily detected in the infrared. 

Even at a resolution of ~5-6", star forming regions 
complexes are sufficiently resolved that their shape have 
to be taken into account using polygonal apertures. 
Visual inspection shows that the irregular shapes are 
mainly due to the overlap of very close star forming re- 
gions. The apertures have been made as small as possible 
to reduce the number of star forming regions inside an in- 
dividual aperture to limit the contamination by nearby 
sources. The goal is to select star forming regions as 
simple as possible, and to avoid the most complex ones 
composed of stellar populations of different ages. To 
achieve this goal, we have defined apertures using poly- 
gons. Each polygon is defined manually to encompass 
as few star forming regions as possible within a single 
aperture as well as to avoid artifacts, background and 
foreground objects. 

Even if some regions are most likely selected in the UV 








Fig. 2. — Combined FUV (blue channel), Ha (green channel) 
and 24 fim (red channel). The top-left image in NGC 5194, the 
top-right one NGC 5236, the bottom-left one NGC 7793 and the 
bottom-right one NGC 5457. 

and in the IR but with different apertures, no attempt 
has been made to make a link between the two. This 
is made difficult by the fact that for a given region the 
apertures can be quite different, cover only a part of the 
region or be split into several apertures. However, a vi- 
sual inspection of the apertures shows that some regions 
have an aperture defined only from either UV or IR even 
if it is visible in both bands. This is due to the fact that 
the apertures for each band are built mainly for the most 
prominent regions. 

The background flux is defined on a per-aperture basis 
using an annulus around the aperture. The radius and 
the width of this annulus is defined for each galaxy de- 
pending on the surface density of star forming regions. 
Local background subtraction has been performed in all 
bands. The background level is then calculated using the 
mode of the pixels enclosed in the background annulus. 
To improve the determination of the galactic background 
level, for some wavelengths we have invalidated emission 
above a given threshold in the background annulus. This 
threshold has been determined manually to eliminate as 
much emission as possible from the nearby sources while 
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TABLE 2 

Images 



Galaxy 


Band 


Telescope 


Date 


Exposure time 

00 


1-cr sensitivity 
(njy arcsec _2 ) a 


NGC 0628 


FUV 


GALEX 


2003-10-07 


1636 


31 




NUV 








28 




U 


Bok 2.3 m 


2004-12-12 


600 


12 




B 


CTIO 1.5m 


2001-10-20 


480 


50 




R 




2001-10-20 


180 


189 




Ha 




2001-10-21 


600 


1.36 x 10" 20 




8 fim 


Spitzer/IRAC 


2004-07-28 


240 


1245 




24 (im 


Spitzer/MIPS 


2005-01-23 


160 


2515 


NGC 2403 


FUV 


GALEX 


2003-12-05 


3317 


22 




NUV 








18 




U 


Bok 2.3 m 


2004-12-09 


600 


35 




B 




2004-12-10 


300 


17 




R 


KPNO 2.1m 


2001-11-16 


180 


74 




Ha 




2001-11-15 


600 


6.70 x 10~ 21 




8 fim 


Spitzer/IRAC 


2004-10-12 


240 


790 




24 (im 


Spitzer/MIPS 


2004-10-12 


160 


1404 


NGC 3031 


FUV 


GALEX 


2006-01-05 


11912 


23 




NUV 




2005-01-12 




17 




U 


Bok 2.3 m 


2004-12-12 


600 


42 




B 




2004-05-10 


300 


48 




R 


KPNO 2.1m 


2001-03-30 


420 


146 




Ha 




2001-03-30 


900 


1.33 x 10^ 20 




8 fim 


Spitzer/IRAC 


2004-05-01 


240 


1219 




24 (im 


Spitzer/MIPS 


2004-10-16 


160 


1597 


NGC 5055 


FUV 


GALEX 


2004-03-31 


1655 


23 




NUV 








23 




U 


Bok 2.3 m 


2005-02-09 


600 


17 




B 






300 


60 




Ha 


KPNO 0.9m 


1999-04-20 


420 






8 fim 


Spitzer/IRAC 


2004-12-15 


240 


1399 




24 (im 


Spitzer/MIPS 


2005-01-26 


160 


1750 


NGC 5194 


FUV 


GALEX 


2003-06-19 


2520 


27 




NUV 








20 




U 


Bok 2.3 m 


2004-06-20 


600 


20 




B 


KPNO 2.1m 


2001-03-28 


720 


37 




R 




2001-03-28 


360 


61 




Ha 




2001-03-28 


900 


5.04 x 10~ 21 




8 fim 


Spitzer/IRAC 


2004-05-22 


240 


1007 




24 (im 


Spitzer/MIPS 


2004-06-22 


160 


1701 


NGC 5236 


FUV 


GALEX 


2003-06-07 


1343 


50 




NUV 








35 




U 


Bok 2.3m 


2005-05-09 


600 


88 




B 






300 


104 




R 




2005-05-10 


300 


125 




Ha 




2005-05-10 


600 


1.37 x 10" 20 




8 /im 


Spitzer/IRAC 


2008-03-05 


240 


1496 




24 (im 


Spitzer/MIPS 


2006-02-18 


160 


2005 


NGC 5457 


FUV 


GALEX 


2003-06-20 


1040 


38 




NUV 








30 




U 


Bok 2.3 m 


2005-05-09 


600 


49 




B 




2005-05-10 


300 


31 




8 /im 


Spitzer/IRAC 


2004-03-08 


85 


1212 




24 (im 


Spitzer/MIPS 


2004-05-10 


200 


1398 


NGC 7793 


FUV 


GALEX 


2003-09-13 


1553 


26 




NUV 








26 




U 


Danish 1.54m 


1997-06-09 


900 


63 




B 


CTIO 1.5m 


2001-10-18 


480 


73 




R 




2001-10-18 


180 


143 




Ha 




2001-10-18 


600 


1.43 x 10" 20 




8 /im 


Spitzer/IRAC 


2004-06-10 


240 


1115 




24 (im 


Spitzer/MIPS 


2004-07-10 


160 


1829 



a In units of W m~ 2 s -1 arcscc~ 2 for Hn images. 
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keeping the genuine background emission. This was done 
for the GALEX bands as well as the Spitzcr IRAC 8.0 /im 
and MIPS 24 ^m bands. This threshold was not applied 
to broadband optical images where young star forming 
regions are much less prominent ; in the U, B and R 
bands the local background is dominated by an evolved 
stellar population. 

The raw fluxes and the errors have been calculated 
using the IRAF procedure POLYPHOT. 

Following the aforementioned criteria, we have selected 
a total of 1146 star forming regions in the 8 galaxies con- 
stituting the sample. In order to increase the significance 
of the results we have subsequently selected only star 
forming regions which have flux uncertainties of at most 
20% in the GALEX FUV, NUV and the Spitzer MIPS 
24 pan. band, leaving us with a total sample of 324 re- 
gions. The number of UV selected and infrared selected 
regions is presented in Table [3l for each galaxy in our 
sample. 

4.2. Aperture correction 

We have performed an aperture correction on the in- 
frared fluxes. As the method routinely used is not tai- 
lored for no n circular aper t ures, we have proceeded as 
described in iBoquien et al.l (|2007l ). We have calculated 
an "equivalent radius" : r = \J S/tt, r being the equiva- 
lent radius and S the surface of the aperture in units of 
pixels. The flux has then been corrected using this equiv- 
alent radius combined with the formulas provided in the 
Sp itzer handbook for exten ded sources for IRAC 9 and 
in lEngelbracht et all (|2007f ) for MIPS. For the 8.0 /an 
(24^m) band the average correction factor was 0.94±0.17 
(resp. 1.70±2.11). No aperture correction has been per- 
formed neither on UV images - as the corrections are 
very small - nor on optical images, thanks to their higher 
angular resolution. 

4.3. Correction of the Galactic extinction 

We have corrected our UV and optical data for 
the effect of fo regrou nd Galactic extinction using the 
ICardelli et al.l (| 1989T ) attenuation curve along with 
the differential extinction E (B — V) provided by the 
NASA/IPAC extragalactic database. 

4.4. Fluxes and main properties of the sample 

The observed fluxes in all bands are presented in Ta- 
bled! The luminosities of the star forming regions span a 
factor -200 in FUV (5.15 x 10 32 to 1.00 x 10 35 W) and a 
factor ~2000 in 24 fim (5.81 x 10 31 to 1.23 x 10 35 W). We 
can see in Figure [3] that those two star formation tracers 
are strongly correlated with a Pearson correlation coeffi- 
cient r = 0.75 between log (L (24)) and log (L (FUV)). 

We observe that especially at higher luminosities, the 
24 /zm luminosity increases faster than the FUV one. 
This is most likely an extinction effect, as can be seen 
in infrared luminous galaxies for instance. The fit of the 
sample gives log L (FUV) = 16.55±0.85+ (0.52 ± 0.03) x 
log L (24). The slope is as expected significantly smaller 
than 1. 

9 Also available at the following address: 
http : //ssc . spitzer . caltech . edu/ irac/calib/extcal/ 




log(L(24)) [W] 

Fig. 3. — Far ultraviolet luminosity versus the 24 (im one. The re- 
gions selected in ultraviolet (resp. mid-infrared) are represented by 
diamonds (resp. by asterisks). The color indicates in which galaxy 
the star forming regions are located: NGC 0628: red ; NGC 2403: 
brown ; NGC 3031: green ; NGC 5055: blue ; NGC 5194: magenta 
; NGC 5236: Orange ; NGC 5457: black ; NGC 7793: cyan. The 
black solid line represents the fit for the whole sample. 

While the correlation between the FUV and the 24 /jm 
star formation tracers is good, the UV bands are im- 
pacted by the dust. The distribution of UV colors is 
presented in Figure [4] 




0.1 0.2 0.3 

log(F,,(FUV)/F,(NUV)) 

Fig. 4. — Distribution of the UV colors as defined by 
log.F\ (FUV) /F\ (NUV). The dashed line histogram represents 
the whole sample. The colored lines represent the histogram for 
each galaxy. The bin width is 0.04 dex. Histograms are slightly 
offset horizontally from each other by a small fraction of their bin 
width for better readability. The colors are the same as defined in 
Figure [3] 

There is a clear distinction in the UV colors distri- 
bution among the galaxies. The bluest one NGC 5457 
peaks 0.2 dex bluer than the reddest one, NGC 5236. 
This may be a consequence of the presence of older UV- 
emitting stars unrelated to the current star formation 
episode or the extinction which affects UV luminosities 
(a range of 0.2 dex corresponds to an extinction range 
Ay ~ 1 for the stellar component assuming a starburst 
attenuation curve). This is supported by the fact that 
very few star forming regions have an UV color that is 
compatible with an ionizing population even though Ha 
is detected in most of them which indicates a redden- 
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TABLE 3 

Number of selected star forming regions 



Galaxy UV selected regions IR selected regions Regions 





Total 


Selected 


Total 


Selected 


Total 


Selected 


NGC 0628 


57 


26 


36 


8 


93 


34 


NGC 2403 


SO 


28 


24 


12 


104 


40 


NGC 3031 


135 


23 


64 


13 


199 


36 


NGC 5055 


51 


7 


40 


5 


91 


12 


NGC 5194 


74 


32 


63 


25 


137 


57 


NGC 5236 


67 


30 


52 


22 


119 


52 


NGC 5457 


143 


50 


121 


34 


267 


84 


NGC 7793 


66 


6 


70 


3 


136 


9 


Total 


673 


202 


473 


122 


1146 


324 



TABLE 4 
Apertures positions and fluxes 





ID 








b 




5* 






FUV C 




NUV C 


U c 


B c 






R c 






8 nm do 




24 fl m d 


Ha" 




R e 










(J2000.0) 


(J2000.0) 


(mJy) 


(mjy) 


(m-Iy) 


(mjy) 




(mjy) 




(mjy) 




(m.Jy) 


10" 


- 17 w m- 


-2 B -l 




NGC 


0628 


U 


01 


01 


36 


28.0 


+ 15 


46 


59 


0.14± 


0.02 


0.13± 


0.03 


0.20± 0.05 


0.23± 





07 


0.17+ 





13 


2.81± 


52 


4.52 + 


0.66 


3.74+ 


0.74 


15.3 


NGC 


0628 


U 


02 


01 


36 


29.1 


+ 15 


48 


21 


0.19± 


o.oi 


0.21± 


0.02 


0.38± 0.04 


0.34+ 





06 


0.37± 





12 


11.96+ 


10 


32.47+ 


0.61 


11.86+ 


0.62 


10.0 


NGC 


0628 


u 


03 


01 


36 


30.3 


+ 15 


■19 


22 


0.29± 


0.05 


0.35± 


0.05 


0.86± 0.10 


1.64± 





11 


1.94± 





30 


7.24+ 1 


)1 


8.67+ 


1.48 


5.49+ 


1.68 


11.9 


NGC 


0628 


u 


04 


01 


30 


36.1 


+ 15 


48 


06 


0.17± 


0.03 


0.22± 


0.03 


0.31± 0.06 


0.41± 





11 


0.36± 





15 


6.30+ 


89 


9.81 + 


1.12 


3.06+ 


0.74 


7.3 


NGC 


0628 


u 


05 


01 


36 


36.2 


+ 15 


47 


26 


0.42 + 


0.07 


0.48± 


0.08 


0.69± 0.18 


0.97± 





39 


0.90± 


0.60 


21.41+ 2 


3 2 


18.28 + 


2.04 


5.95 + 


1.67 


12.3 


NGC 


0628 


u 


06 


01 


36 


36.4 


+ 15 


50 


17 


0.55± 


0.10 


0.57± 


0.10 


0.95± 0.18 


1.03± 





31 


0.96± 





53 


13.90+ 2 


61 


24.61 + 


4.09 


22.13 + 


2.90 


12. 8 


NGC 


06 28 


u 


07 


01 


36 


37.1 


+ 15 


48 


09 


0.35± 


0.05 


0.34± 


0.05 


0.55± 0.12 


0.70± 





2 1 


0.74± 





2!l 


26.09+ 3 


14 


45.01 + 


2.90 


13.02 + 


1.38 


8.1 


NGC 


0628 


u 


08 


01 


36 


37.8 


+ 15 


45 


01 


0.55± 


0.06 


0.66± 


0.07 


0.85± 0.12 


0.95± 





22 


0.82± 





30 


21.96+ 1 


70 


70.66 + 


2.14 


17.26+ 


1.56 


16.2 


NGC 


0628 


u 


09 


01 


30 


38.8 


+ 15 


44 


12 


1.97± 


0.18 


2.14+ 


0.19 


2.60± 0.30 


3.02 + 





58 


2.34± 





81 


26.32+ 2 


94 


76.63 + 


2.64 


37.59+ 


4.51 


8.5 


NGC 


0628 


u 


10 


01 


36 


38.8 


+ 15 


49 


13 


0.69+ 


0.04 


0.83± 


0.04 


1.05± 0.08 


1.50± 





13 


1.16± 





29 


7.95+ 


82 


13.19 + 


0.83 


6.45 + 


1.18 


17.1 



NOTE. Units of right ascension an 

o infrared bands. Galactic extinction 
ournal. A portion is shown here for ; 

Format: names of the galaxy-band i 



; been perform 
ding its form ; 



nd units of dcclinatior. 
on UV and optical flu 



• degrees, arcminutcs, and arcscconds. Aperture correction has only been applied 
Tablc[^is published in its entirety in the electronic edition of the Astrophysical 



U (respectively I) indicates that the apcrtu 



defined from the FUV (resp. 24 fim) : 



1 Posith 

Stellar and 
1 Infrared flu 



of the 



ed ga 



enter 
flu 



ected for the Gala 



vo been corrected for aperture effects. 

. the 8 [ira has been subtracted using the rclati 



: foreground extinction using the Cardclli et a.l. I 1989) law. 

n presented in |rTclou ct al.| [2664) : F, y (8 PAH) = F v (8) - 0.232 X F u (3.6). 



ing unrelat ed to the age o f the current episode. Indeed, 
assuming a iKroupal (| 200 If) IMF (Initial Mass Function), 
and a solar metallicity, a simple stellar population will 
still have EW(Ha)~lA after 20 x 10 6 years. 

Among the 324 selected star forming regions, we have 
suitable Ha imaging information for 228 and we measure 
their fluxes. The lack of Ha observations for NGC 5055 
(we do not have the corresponding suitable R band im- 
age to perform continuum subtraction making the Ha 
fluxes too unreliable to be used here) and NGC 5457 
is responsible for this difference. Among those regions, 
213 (respectively 184) are detected with S/N > 3 (resp. 
S/N > 5). It shows that the star forming regions are 
young and have recently formed stars. Indeed, the Ha 
equivalent width is a good tracer of the age of a single 
star forming episode for the first ~ 10 — 15 x 10 6 years. 
We will see in section 16.21 that the U — B color is also a 
good tracer of the age. 



4.5. Estimation of the extinction 

In order to estimate attenuation values for our regions, 
we ha ve made use o f the 24 /zm and the Ha images. 
Indeed iCalzetti et al.l (|2007l ) showed that there is a very 
good correlation between the extinction corrected Paa 
and a combination of the 24 fim. and dust obscuratcd Ha 
leading to an estimate of the intrinsic Ha luminosity: 



L (Ha) 



ved 



0.031 xL (24). (1) 



This permits to derive the extinction of the gas: 



j^gas 



3.06 x log 



1 



0.031 x £(24) 
MHa) 0l 



bserved 



(2) 



As the extinction of the stars is smaller, we assume here 
that A v = 0.44 x A 9 y S (jCalzetti et al.ll2000h . We obtain 
an extinction ranging from Ay ~ 0.06 to Ay — 1.28 with 
a mean {Ay} = 0.44 ± 0.28. 

Note that as this method to estimate the extinction is 
dependent on the Ha, NGC 5055 and NGC 5457 fluxes 
cannot be corrected for the extinction and arc therefore 
excluded from any subsequent analysis involving extinc- 
tion correction. 

5. MODELING 
5.1. Spectral synthesis 

To model the SED of individual star forming regions, 
we have use d the Pegase ii evolutionary sp ectral syn- 
thesis code (jFioc k, Rocca-Vohnerangel 119971 ). We have 
assumed an instantaneous starburst having an age rang- 
ing from 1 x 10 6 years to 300 x 10 6 ye a rs wit h a constant 
metallicity of Z = 0.02 and a IKroupal (|2001[ ) initial mass 
function from 0.1 M Q to 120 M©. The impact of assum- 
ing solar metallicity in the models on our results will bc 
examined in section f6. 21 

Fluxes and colors are determined convolving the gen- 
erated SED - taking into account the nebular emission 
lines as well as the nebular continuum which is important 
for young regions (|Anders fc Fritze-v. Alvenslebenll2003T) 
- with the spectral response function of the filters for the 
FUV, NUV, U and B bands. The infrared luminosity is 
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obtained calculating the absorbed flux by the dust. See 
the section hereafter for details about the dust models. 

5.2. Dust models 

We use different dust models in order to probe the 
effect of dust- stars geometry: 1) the starburst atten- 
uation curve (|Calzetti et al.l [2000h : classically used in 
deep surveys with the IRX-/3 relation, 2) the foreground 
dust screen Small M agellanic Cloud attenuation curve 
([Gordon et al.l[2003h . 3) the mixed stars and dust Small 
Magellanic Cloud attenuation curve using the relation: 

WrvedA ntrmsic = 0-921 X (l - 1(T°- 4 ^)) /A (A) 

adapted from lNatta fc Panagial (|1984f l. 

The dust extinction is applied to the unextinguished 
SED with Ay ranging from to 3 by steps of 0.1. In 
this paper, conventionally, Ay designates the extinction 
of the stars. The extinction of the gas is noted A 9 y h . The 
nebular emission accounts for a small fraction of the total 
flux in the broad bands, and the same extinction law has 
been applied to the stellar and the nebular emissions in 
these cases. 

A comparison with the models generated by 
iCalzetti et al.1 (|2005l ) using Starburst99 shows that the 
results are nearly identical. 

6. THE AGE AS THE "SECOND PARAMETER"? 

In order to probe the driving parameters for the spread 
of the normal star-forming galaxies in the IRX-/3 dia- 
gram, we need to derive the following quantities for our 
galaxies: the ultraviolet slope, the total infrared lumi- 
nosity (section 16. ip , which will be used to derive the 
FIR/FUV ratio, and estimators for the age of the stellar 
populations. In the present work, we use two estimators: 
the U — B color and the equivalent width of Hq emission 
line (section [6. 2|) . 

6.1. Estimation of the bolometric infrared luminosity 

In order to determine the total infrared emission of 
the selected star forming region we use the Spitzer IRAC 
8.0 /jm and MIPS 24 (im emission. While the exploita- 
tion of the MIPS 70 /im and 160 fim bands would yield 
better results, their resolution (17" and 38") is too low 
to separate individual star forming regions. For this 
study we estimate th e total infrared emission using the 
ICalzetti et"all (120051 ) relation that was derived in a pre- 
liminary study on NGC 5194: 

log L (IR) = log [L (24)]+0.908+0.7931og[£ 1/ (8) /L v (24)] . 

The scatte r around the relatio n within a galaxy is 
about 40% (|Calzetti et alJ 120051 ). which is the accu- 
racy level we can expect for the determination of 
the total infrared flux of ind iv idual star forming re- 
gions. IPerez-Gonzalez et al.l (|2006l ) also derived a 
sim ilar relation. Pract ically, th e difference between 
the | Calzetti et all (|2005l ) and the IPerez-Gonzalez et al.1 
(l200l ~ IRX estimates for the selected star forming re- 
gions is —0.07 ± 0.03 dex which is only a few percent of 
the total range and does not alter the results presented 
in this study. 

6.2. Stellar population age estimators 



To test if the age is the "second parameter" we are 
looking for, an accurate determination of the stellar pop- 
ulation ages is needed. Given the large number of star 
forming regions involved, a direct spectral observation of 
the age sensitive 4000 A break, as measured by the -D4000 
parameter are not available for the majority of regions in 
our sample. We have decided to exploit the U — B color 
as a proxy to the 4000 A break. Indeed the U band is cen- 
tered around 360 nm while the B band is centered around 
440 nm, on each side of the bre ak. A simple Pegase 
II model using a IKroupal (|2001l) IMF from 0.1 M Q to 
120 Mq and a metallicity of Z=0.02 shows that for an 
instantaneous burst, the U — B color ranges from —1.18 
after 1 x 10 6 years to 0.04 after 300 x 10 6 years (Figured] 
left). 

While the -D4000 parameter is essentially extinction 
free thanks to the narrow wavelength range probed, ex- 
tinction is strong at short wavelengths, inducing an im- 
portant uncertainty regarding the broadband U and B 
fluxes. However this problem is mitigated when using 
the U — B color. Indeed, when averaged on a number 
of dust geometries and extinction curves, an extinction 
of Ay = 1 yields only E (U - B) ~ 0.1 x Ay mag. We 
will see in section 16.31 that the stellar extinction - as- 
suming A^ ars = 0.44 x A 9 y s - ranges from Ay = 0.06 
and Ay = 1.28, with a mean extinction of 0.44 magni- 
tude, which yields an average uncertainty of 0.04 magni- 
tude (0.1 magnitude if we assume a starburst attenuation 
law). 

Young star forming regions are the main source of 
Ha emission in star forming galaxies. As the popula- 
tion ages, the Ha drops dramatically as stars emitting 
ionizing photons have a life expectancy typically under 
20 x 10 6 years. A standard way to estimate the age is 
the equivalent width of the Ha emission line, which de- 
creases as the most massive stars die. The exact esti- 
mation would necessitate spectroscopic observations. In- 
deed, even if the R band emission is not constituted only 
of the continuum - it is "polluted" by the Ha and [Nil] 
lines and the nebular emission in general - it is the best 
possible approximation with the available data. Another 
advantage is to limit the differential extinction as the Ha 
line and the R band central wavelength are close. 

A comparison of the evolution of U — B and Ha equiva- 
lent width age tracers against time is plotted in Figure [5] 
left and they are compared with observations in Figure [5] 
right. 

We see in Figure [5] right that for all star forming re- 
gions, the Ha equivalent width is reasonably high, rang- 
ing typically from ~ 15 A to ~ 600 A. This indicates 
a young age, typically less than 10 x 10 6 years. The 
U — B color ranges typically from —1.3 to —0.1. There 
is a clear correlation between the Ha equivalent width 
and the U — B color with a Pearson correlation coeffi- 
cient r = —0.52. NGC 5194 star forming regions are 
significantly offset from the mean correlation, having a 
much lower Ha equivalent width for their U — B color (or 
having a much bluer U — B color for their Ha eq uivalent 
width ). The instantaneous burst model with a iKroupU 
(|200lD IMF and a metallicity of Z = 0.02 can reproduce 
most star forming regions that have a U — B color redder 
than —0.9 assuming an extinction such as E (U — B) = 
0.24 x Ay mag and EW (Ha) oc iq-oaxo.68xA v with 
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1 10 100 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 

t [Myr] U-B 



Fig. 5. — Left: U — B co l or (bl ue line) and the Ha equivalent width (red line) versus the time elapsed since an extinction- free instantaneous 
starburst with a IKroupal (120011 ) IMF and a metallicity Z = 0.004 (dashed line), Z = 0.02 (solid line) and Z = 0.01 (dotted line). Note 
that practically measurements of equivalent widths smaller than a few A are uncertain. Right: equivalent width of the Ho line defined as 
L (Ha) /L\ (R) in units of A. The colors are the same as defined in Figure [3] The solid colored line represents the track of an extinction- 
free instantaneous starburst with a[Kroupa (200l|) IMF and a metallicity Z = 0.02. The age ranges from 10 6 years to 20 X 10 6 years by 
increments of 10 6 years. It is indicated by the colorbar. The black solid lines departing from the colored solid line represent the effect of 
the extinction, assuming that E (U — B) = 0.24 X Ay mag and EW (Ha) oc 10~ 0-4x 0.68 x .Ay ma g which takes into account the differential 
extinction between the gas and the R band emitting population. Each intersecting line represents an increment of the extinction by 0.1 
magnitude. 



0.04 < A v < 1.28. This shows that while the U-B color 
is only weakly sensitive to the extinction, the inclusion 
of extinction effects is nevertheless needed to reproduce 
most of the observations. However for the bluest U — B 
star forming regions, the model significantly overpredicts 
the observations of the Ha equivalent width. One pos- 
sible explanation is an insufficient background removal 
in the R band image within each region, artificially de- 
creasing the EW (Ha) values for blue U — B colors. An 
alternative explanation is a strong contamination of the 
R band by the Ha and [Nil] lines which would lead to 
significantly overestimate the continuum in the R band. 
However, the mean EW(Ha) is 195 ± 193 A and the typ- 
ical width of a R band filter is ~ 1500 A. This means the 
contamination is at most about 25%, insufficient to re- 
produce the observations. In section [6731 we will see that 
the characteristics of most of the star forming can be re- 
produced by a simple model in the IRX-/3 diagram. An 
important fraction of them however should not present 
any significant Ha emission - 73% (rcsp. 60%) should 
be older than 12 x 10 6 years (resp. 20 x 10 6 years) - 
if the underlying hypotheses (the star formation history 
can be modeled as an instantaneous burst and the trend 
observed as a function of the perpendicular distance is 
an age effect) are correct. 

A parameter that may affect the age estimate is the 
metallicity. Indeed, higher metallicity stars have a sig- 
nificantly redder U — B color due to the presence of a 
high number of metallic absorption lines in their atmo- 
sphere. An illustration of this effect on the U — B and 
the ultraviolet colors is presented in Figure [6] 

While the metallicity has an influence, it is small as the 
oxygen abundances are typically in the relatively narrow 
range 12 + logO/H =8.6-9.2, except NGC 5457 which 
outermost r egions have a metalli city down to 8.2. As 
a reference, lAsplund et alj (|2005f ) determined the solar 
metallicity as Z = 0.0122 or 12 + logO/H = 8.66. 

Another uncertainty is due to the synthesis of the U- 




1 10 100 

t [Myr] 

Fig. 6. — Effect of the metallicity on the U — B and ultra- 
violet colors. The solid lines represent the U — B colors while 
the dashed ones represent the log (Fx (FUV) /F\ (NUV)) colors. 
The red lines give the difference between Z = 0.02 and Z = 0.1 
while the the blue ones give the difference between Z = 0.02 
and Z = 0.004. For an ionizing population (17 — ^)z— o l — 
(U-B) Z=Q0M (resp. log (F x (FUV) /F x (NUV)) Z=0 1 - 
log (F x (FUV) /F x (NUV)) z=omA ) reaches a maximum of 0.26 
(resp. 0.16). 

band luminosities. Indeed, even in simple U — B, B — V 
color-color diagrams, the model tracks are not consistent 
with integrated photometry of dwarf galaxies (|Ledl2006l . 
for more details). 

6.3. IRX-(3 diagram 

In Figure [7] we present the IRX-/3 diagram of the se- 
lected star forming regions compared with the model pre- 
sented in section [5] 

Several results have to be noted. We see that most of 
the star forming regions are in general agreement with 
the models though a number of regions have a conspicu- 
ously blue UV color or a particularly low IRX. We also 



10 




Fig. 7. — IRX-/3 diagram of individual star forming regions. The regions represented by a diamond (resp. star) are FUV (resp. 24 /tun) 
selected. The color of the symbols indicates the actual U — B color (i.e. the age) following the scale defined by the color bar on the right 
corner. The error bars are plotted in light gray. The grid is the model described in section [5] using a starburst attenuation curve (left) and 
the SMC extinction curve (right). The color of each vertex corresponds to the expected U — B color. The axis predominantly horizontal 
(resp. vertical) samples the age: 1 X 10 6 , 6 X 10 6 , 12 X 10 6 , 20 X 10 6 , 30 X 10 6 , 40 X 10 6 , 50 X 10 6 , 75 X 10 6 , 100 X 10 6 , 150 X 10 6 , 
200 X 10 6 , 250 X 10 6 and 300 X 10 6 years (resp. extinction from Ay = 0.1 to Ay = 3.0, the first horizontal line representing an extinction 
of 0.1 magnitude.). The solid lines represent an ionizing population while the dashed lines represent a non-ionizing population. The model 
representing the origin of the perpendicular distance (d = 0) corresponds to an age of 6 X 10 6 year s. We see on the grid t hat the age has a 
much stronger effect on the U — B color than the extinction. The black squares represent SINGS (Kcnnicutt ct al. 2003) galaxies and the 
black triangles LVL ones (Dale et al., 2009). 



see that the selected star forming regions describe a clear 
envelope mostly circumscribed in its upper part by the 
modeled starburst IRX-/3 relation of a very young star- 
burst. We retrieve here the well known result that quies- 
cent star forming galaxies lay under the starburst IRX-/3 
curve. 

A visual inspection shows an expected qualitative 
trend to have younger regions where the UV color is the 
bluest while on the right part there are predominantly 
older regions which have a redder UV color. Interestingly 
the younger regions have a trend closer to the modeled 
IRX-/3 starburst relation. One explanation could be that 
it marks the locus of both constant star formation and 
young stellar regions. Indeed, in addition to their young 
age, the dust-stars geometry is likely to be similar to 
that of starburst galaxies (star forming regions inside a 
dust shell). Another explanation could be that those re- 
gions are still relativel y enshrouded in their birth cloud. 
However, as shown bv lPrescott et al.l (|2007|) . only 3% of 
regions are completely enshrouded in dust. 

6.4. Age sensitive tracers as a function of the 
perpendicular distance 

Here we quantify whether age may be the parameter 
accounting for the deviations of our regions from the star- 
burst attenuation curve in the IR-/3 diagram. 

A standard tool to evaluate the distance of a star form- 
ing region (or a galaxy as a whole) from the standard 
starburst attenuation l aw is the perpend icular distance 
in the IRX-/3 diagram (K ong et al.l 120041 , for instance). 
The distance is simply defined as the euclidean distance: 

d = \J (x - x E l) 2 + {y- Vel) 2 , (4) 

where x (resp. xel) and y (resp. jael) are the /3 and 
IRX values of the individual star forming regions (resp. 
of the extinction law locus that is the closest to the ob- 
served value) and d is the dimensionless distance in the 



IRX-/3 plot. The starburst and SMC IRX-/3 relations are 
obtained from the model presented above. The SMC ex- 
tinction law will be studied more in detail in section [71 
We assume an age of 6 x 10 6 years which gives an intrinsic 
log F x (FUV) /F x (NUV) = 0.37. The IRX-/3 relation is 
obtained varying th e extinction applie d to this popula- 
tion. Following the iKong et all (|2004fl convention, star 
forming regions which have a higher (resp. lower) than 
expected IRX given their UV slope will have d > (resp. 
d < 0). 

In Figure [5] we plot the equivalent width as the ratio 
between the Ha flux and the R band flux density cor- 
rected for the extinction using the formula in section [6~2l 
against the perpendicular distance. 

As expected the mean Ha equivalent width increases 
as a function of the perpendicular distance with a slope 
of 385 and a large la scatter around the linear fit: 150A. 
The perpendicular distance and the Ha equivalent width 
are weakly correlated with a Pearson correlation coeffi- 
cient r = 0.18. The most surprising result is that even at 
large negative perpendicular distances, a significant Ha 
emission is detected. This result shows that the age as 
measured by the Ha equivalent width is not the "second 
parameter" we are seeking. From the point of view of 
star formation, it could mean that the assumption of an 
instantaneous starburst is not valid. We will explore if 
variations of the star formation history can explain the 
observations in section f6. 51 

In Figure [9] the extinction corrected U — B - assuming 
a starburst attenuation law: E (U — B) = 0.24 x Ay mag 
- is plotted as a function of the perpendicular distance. 

While a trend of the U — B color as a function of the 
perpendicular distance could be expected, it is partic- 
ularly weak, with a slope of only —0.62 ± 0.20. The 
perpendicular distance and the U — B color are corre- 
lated with a Pearson correlation coefficient r = —0.17. 
The weakness of the trend combined with the important 
scatter hints that another physical phenomenon than the 
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Fig. 8. — Differential extinction corrected Ho equivalent width , as d efined by the ratio of the Ha over R band flux density, as a function of 
the perpendicular distance (as defined in equation [4] see section [6,41 1 between the starburst (left) and SMC (right) laws and the individual 
star-forming regions in the IRX-/3 diagram. The purple line is the fit of the equivalent width as a function of the perpendicular distance: 
EW (Ha) = 264 ± 13 + (385 ± 144) x d (starburst) and EW (Ha) = 234 ± 10 + (410 ± 157) x d (SMC). 
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Fig. 9. — Extinction corrected U — B color as a function of the perpendicular distance (as defined in equation|4] see section [6.41 1 between 
the starburst (left) and SMC (right) laws and the individual star-forming regions in the IRX-/3 diagram. The purple line is the linear fit 
of the extinction corrected U — B color as a function of the perpendicular distance: U — B = —0.85 ± 0.02 — (0.62 ± 0.20) X d (starburst) 
and U — B = -0.80 ± 0.01 - (0.74 ± 0.22) x d (SMC). 



age plays a much more important role. 

We have shown here that the age, which was a prime 
candidate to explain the deviation from the starburst 
attenuation law in the IRX-/3 diagram does not seem 
to play a major role. Indeed, using two different age 
tracers, the U—B color and the Ha line equivalent width, 
the large dispersion around the relation shows that the 
influence of the age is at most minimal under the current 
assumptions and for the timescales probed by those two 
tracers. 

6.5. Influence of the star formation history and dilution 
by an older stellar population 

We have shown that an instantaneous starburst cannot 
reproduce simultaneously the large perpendicular dis- 
tances and a non negligible Ha emission. In particular, 
there are a number of regions in our sample that are too 
red in UV colors. We now consider the possibility that 
the stellar populations in those star forming regions are 
more complex than a simple instantaneous burst descrip- 
tion. We thus introduce a model of the stellar population 



of a galactic disk that considers an exponentially decreas- 
ing SFR of the form: 

SFR (t) = a x exp (-i/r) , (5) 

where a is an arbitrary normalization constant in units of 
M0.yr _1 , t is the time in years and r is the time constant 
in years. 

First of all we want to see if an exponentially de- 
creasing star formation rate is able to reproduce the 
observations. To take into account different cases of 
star formation history we test 3 models which differ 
by the ratio SFR (12 x 10 9 ) /SFR{0). SFR(0) is the 

SFR at t = year and SFR (12 x 10 9 ) is the SFR at 
t = 12 x 10 9 years, both from equation [5] The values of 
this ratio are 1/10, 1/100 and 1/1000 (which correspond 
to a birthrate parameter of 2.6 x 10 _1 , 4.7 x 10~ 2 and 
6.9 x 10~ 3 and a timescale of 5.2 x 10 9 years, 2.6 x 10 9 
years and 1.73 x 10 9 years respectively). It turns out that 
this star formation history can be readily excluded as it 
cannot reproduce at all the observations. The reason is 
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that even for SFR (12 x 10 9 ) /SFR(0) = 1/1000, the 
current star formation rate is high enough to overwhelm 
ultraviolet emission from older stellar populations. 

A more realistic model for a galactic star forming re- 
gion is a small burst on top of the galactic stellar popu- 
lation. To model the galactic stellar population we use 
an exponentially decreasing star formation rate as de- 
scribed in equation [5] for < t < 11.95 x 10 9 years, no 
star formation for 11.95 x 10 9 < t < 12 x 10 9 years and 
an instantaneous burst at t = 12 x 10 9 years. The ratio- 
nale is that in the area of the star forming region, due to 
the diffusion of the stars over time, the stellar population 
is made of stars born in different regions of the galaxy 
subsequently mixed. However stars younger than a few 
10 x 10 6 years are too young to have migrated long dis- 
tances. In addition we assume that in the last 50 x 10 6 
years, no star formation has taken place at the location 
of the current star forming region and no significant mi- 
gration from nearby star forming regions has occurred. 
The mass of the instantaneous starburst at t = 12 x 10 9 
years is defined as M stQrs /300 following iCalzetti et al.1 
were M star s is the mass of the stars constituting 
the background galactic population. The model is plot- 
ted in Figure [10] along with the observations in an IRX-/3 
diagram. 
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Fig. 10. — Same as Figure [7] with a model constituted of an 
exponentially decreasing burst during 11.95 X 10 12 years and an 
instantaneous starburst of a mass of M st ars/300 at t = 12 x 10 9 . 

While this model reproduces a much larger range of 
ultraviolet colors as a whole, it still cannot explain the 
Ha emission of such regions. Indeed, initially the in- 
stantaneous starburst emission dominates over the older 
population until it fades away after a few tens millions 
years, after which the older population dominates and is 
responsible of the redder ultraviolet colors compared to 
the model presented in section 

While a continuous star formation does not account 
for the observations, we can take into account a more 
stochastic nature of the star formation history by con- 
sidering a single, still UV emitting, older burst. Indeed, 
the presence of a significant old stellar population can 
have an influence on the UV color. Individual star form- 
ing regions with local backgrounds subtracted are less 
affected by this effect than integrated galaxies. Some 
of the observed regions are redder than can what can be 



reproduced by an ionizing population and a SMC extinc- 
tion law (see Figure [7] and section [3 the SMC extinction 
law is already extreme in that it has a very steep slope 
in the UV, this means that for a given IRX, it reddens 
more the UV color than any other extinction law consid- 
ered in this study). That there are some regions which 
have an even redder color hints about the possibility of 
a dilution by older stellar populations. In order to test 
this, we have made a simple model assuming 2 instan- 
taneous bursts, the one we are currently observing and 
another one which is 300 x 10 6 years old. The star forma- 
tion history is undoubtedly more complex but cannot be 
precisely retrieved with the current observations. How- 
ever, NUV emitting populations have an life expectancy 
of about 10 9 years. We take here a 300 x 10 6 years old 
instantaneous burst as an approximation of the real star 
formation history. The mass ratio of the more evolved 
stellar population to the younger one varies from to 
300. We have selected all star forming regions that have 
a redder UV color than a 20 x 10 6 years old starburst 
extincted with the SMC extinction law. We then have 
calculated what dilution would be the necessary to red- 
den the UV color, if those regions were actually 20 x 10 6 
years old to account for the fact that they still emit in 
Ha. There is a total of 67 star forming regions which 
have a UV color redder than a SMC-extincted 20 x 10 6 
years old burst. The model shows that the mean mass 
ratio is 7.8 ±7.1 with a median of 5. This is much smaller 
than the ratio for integrated galaxies which reach up to 
few hundreds. This corresponds to a mean flux ratio of 
0.15 in FUV and 0.29 in NUV. Note however that 64% 
of those star forming regions have flux uncertainties such 
that the observations could be reproduced without need- 
ing an older stellar population. 

The dilution by an old stellar population has an influ- 
ence not only on the UV color but also simultaneously 
on the U — B color as we can see in Figure [11] 
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Fig. 11. — Influence of old stellar population on a 10 6 and 20xl0 6 
years old, Ha emitting, burst on the UV and U — B colors. The 
straight solid lines represents the locus of a SMC-extincted 10 6 
(lower line) and 20 X 10 6 years old (upper line) bursts from Ay = 
0.1 to Ay = 1.0. The color, solid, dotted and dashed lines trace 
the evolution of the colors as a function of the mass ratio for an old 
population of 300 X 10 6 , 500 X 10 6 and 10 9 years respectively. The 
successive filled squares (for the 20 X 10 6 years old population) and 
the triangles (for the 10 6 years old population) indicate the mass 
ratio: 10, 25, 50, 100, 150, 200, 300. 
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Wc see that even a relatively small mass ratio can have 
a significant influence on the colors. The influence fades 
away as the older population gets older. Indeed, the 
intrinsic UV luminosity drops quickly with age. Hence, 
even if the older population gets redder, a higher mass 
ratio is needed to offset its smaller UV luminosity. On 
the other hand, the life expectancy of U and B emitting 
stars is much longer than those which predominate the 
UV luminosity. That is why an older population has 
significantly more influence on the U — B color than on 
the UV color. 

7. EXTINCTION LAWS 

We will explore here the effect of varying our assump- 
tions on the dust attenuation law and dust geometry. 
We add to the starbur st attenuation curve the case of an 
SMC extinction curve (jGordon et al.ll2003l ) for both fore- 
ground and mixed dust geometries. We will not c o nsider 
curves like the Milky Way one, as iCalzetti et all ([20051 ) 
showed that this curve cannot reproduce the observed 
trend of the data in the IRX-/3 diagram, for any dust 
geometry. Conversely, the SMC extinction law, coupled 
with a foreground screen seems effective at reproducing 
part of the observed trend (Figure [7|) . 

We observe that most of the star forming regions can 
be reproduced by one or more models but that no model 
can reproduce all observations. Because the SMC extinc- 
tion law has a sharper rise compared to the starburst 
attenuation law, the UV color is more sensitive to the 
extinction in the former case than in the latter one, thus 
producing redder UV colors for equal FIR/UV ratios. 

While the model using the SMC extinction law fails to 
reproduce the youngest regions that could be explained 
with a starburst attenuation curve it can however repro- 
duce a number of UV red regions with a younger inferred 
age. The question of the reproducibility of the Ha emis- 
sion remains in fact open as most of the selected star 
forming regions are detected in Ha. In order to address 
this question we plot the UV colors that can be repro- 
duced by Ha emitting populations as a function of the 
extinction in Figure 1121 For each model we plot two 
lines, one for a 10 6 year old starburst and one for an age 
of 12 x 10 6 years (the oldest age that can still produce a 
reasonably high number of ionizing photons) . 

We observe that as expected, a single extinction law 
cannot reproduce all of the observations. For instance, 
for a given extinction, in the case of the starburst atten- 
uation law, the ultraviolet color is too blue. Regarding 
the SMC extinction law, the behavior is very different 
depending on the assumption made on the stellar extinc- 
tion. If we assume Ay = , only the least extin- 
guished star forming regions can be reproduced. Dimin- 
ishing the extinction of the stars relative to the extinction 
of the gas, more star forming regions can be reproduced 
by the model, but only a fraction of the total number at 
a time (i.e. by each model). 

So far we have always inspected the evolution of age 
tracers assuming the starburst attenuation law. As 
stated above, if the actual extinction law in quiescent 
star forming galaxies is different, it can emulate an age 
effect. As we have already specified, star forming regions 
in the IRX-/3 diagram are bracketed by the starburst at- 
tenuation law and foreground SMC extinction law. 

In Figure [8] right, wc plot the Ha equivalent width 
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Fig. 12. — Ultraviolet color as a function of the extinction in 
magnitudes. Star forming regions are plotted with symbols and 
colors as denned in Figure [7] The starburst attenuation law is 
plotted with a solid line, SMC with a dashed line, SMC(A V /0.44) 
— this is equivalent to Ay = Ay'" , assuming that the relation 
Ay = 0.44 X Ay^ s is also valid in the case of a SMC attenuation 
law - with a dotted line, and a mixed SMC law with a dash triple- 
dot line. The blue lines are for a 10 6 years old burst and the red 
ones for a 12 X 10 6 years old burst. 

versus the perpendicular distance from the SMC IRX-/3 
law. The trend between the Ha equivalent width and the 
perpendicular distance is even weaker than in the case of 
a starburst attenuation curve: the slope is similar (410 
versus 385) as well as the Pearson correlation coefficient 
(r = 0.13 in the SMC case and r = 0.18 in the starburst 
one). 

In Figure [9] right, we plot the U — B age tracer ver- 
sus the perpendicular distance from the SMC IRX-/3 law. 
The trend between the U — B and the perpendicular dis- 
tance is approximately as weak as in the case of a star- 
burst attenuation curve: the slope is similar (—0.74 ver- 
sus —0.62) as well as the Pearson correlation coefficient 
(r = -0.22 versus r = -0.17). 

The standard IRX-/? relation is based on the starburst 
attenuation law. However the physical conditions for 
such a law to be valid are not necessarily always met in 
the case of star forming regions. In this section we have 
investigated different extinction laws. We have shown 
that no single law is able to reproduce the observations. 
Conversely, a range of extinction laws between the star- 
burst attenuation law and the foreground SMC extinc- 
tion law can explain the obser v ations , in line with what 
was found bv lBurgarella et all (|2005f ). 

8. METALLICITY 

The metallicity can significantly affect the ultraviolet 
and U — B colors as we have seen in section 16.21 In 
Figure [13] wc plot 3 models (using the same parameters 
as the ones described in section 1572"]) which differ only by 
the metallicity: Z = 0.1, Z = 0.02 and Z = 0.004. 

The most metal poor model spans a small UV color 
range, slightly more than 0.2 dex. The most metal rich 
model conversely has a much faster evolving ultraviolet 
color. While this may explain the presence of Ha emis- 
sion in a number of regions it nevertheless still cannot 
account for the most UV-red regions, assuming a star- 
burst attenuation law. 

As we see, the metallicity has an important effect on 
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Fig. 13. — Same as Figure [7] The blue grid represents the model 
with Z = 0.004, the green one with Z = 0.02 and the red one with 
Z = 0.1. 

the location of a star forming region on the IRX-/3 di- 
agra m. Another effect o f the mctallicity is the extinc- 
tion (jCalzetti et al.l[T994h . Indeed, it is well known that 
star forming regions with a high er metallicity tend to 
be more extinguished on average (IHeckman et alj 119981 ; 
iBoissier et "all 12004 iMunoz-Mateos et al.ll2009[ ). An ac- 
curate determination of the metallicity of 324 star form- 
ing regions in 8 different galaxies would necessitate long 
spectroscopic observations that are not available. How- 
ever, we also know that there is a metallicity gradient in 
spiral disks. While this does not provide the exact metal- 
licity of a star forming region from its deprojected galac- 
toccntric distance, it statistically provides a precious in- 
formation. In order to evaluate the mctallicity, wc use 
the metallicity gradients provided by Moustakas et al. 
(2009, in preparation) and additional data provided in 
Table [T] The position angles have been calculated from 
optical images. The corresponding IRX-/3 diagram is pre- 
sented in Figure 
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Fig. 14. — IRX-/3 diagram. The color of the star forming regions 
indicate their metallicity. The two curves are IRX-/3 relations for 
the starburst attenuation law (blue) and the foreground SMC (red) 
extinction law. The age of the instantaneous burst is 6 X 10 6 years 
old. 

It appears there is indeed a strong correlation between 
the metallicity and the extinction. The low metallic- 



ity regions have a blue ultraviolet color and a low IRX 
while regions with the highest metallicity have the red- 
dest ultraviolet color and a high IRX. This is consis- 
tent with the models presented in Figure [13] where for a 
given age, a star forming region with a higher mctallic- 
ity would have a redder ultraviolet color. The Pearson 
correlation coefficient between the metallicity and IRX 
is r = 0.80. The relati on is tighter than the one found 
bv lCortese et al.1 (|2006i ) for optically selected, integrated 
normal star-forming galaxies, they derived a Spearman 
correlation coefficient r s = 0.59. 

In order test quantitatively the dependence of the ex- 
tinction on the metallicity, we plot in Figure [T5l the IRX 
value of each star forming region versus the oxygen abun- 
dance 12 + log O/H. 
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Fig. 15. — IRX versus the oxygen abundance determined using 
the galactic mctallicity gradient. 



The Pearson correlation coefficient between the metal- 
licity and IRX is r = 0.80. 

In this section we have examined whether the metal- 
licity can explain the deviation from the IRX-/3 relation. 
We have shown that the metallicity however is an im- 
portant parameter to determine the location of the star 
forming regions in the IRX-/3 diagram. The more metal- 
lic they are, the higher IRX and the redder (3 they have, 
closely following IRX-/3 relations but showing no obvious 
trend with the perpendicular distance. 

9. DISCUSSION 

Understanding the deviation from the starburst atten- 
uation law in the IRX-/3 diagram is essential in our un- 
derstanding of star forming galaxies as they provide im- 
portant insights into the attenuation and reprocessing of 
starlight. We have shown that the age of the star forming 
regions - a debated parameter - has very little influence 
on the deviation, which actually necessitates a range of 
geometries and extinction laws. In addition, the influ- 
ence of the star formation history and in particular any 
dilution by an older population should not be neglected. 
If these parameters are not taken into account properly, 
the error on the determination star formation rates can 
be severely impacted. 

9.1. Age of the star forming regions 
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We have shown in section ROl that there is a weak corre- 
lation between the perpendicular distance from the IRX- 
P relation and the age as traced by the Ha equivalent 
width and the U — B color. In all cases, the trends with 
the age are too weak relative to the scatter in the data to 
support the notion of aging stellar populations as an im- 
portant "second parameter" , at least for the timescales 
probed by the U — B color and the Ha equivalent width. 
We have seen in Figures [8] and [9] that neither the U — B 
color nor the Ha equivalent width correlate strongly with 
the perpendicular distance. The absence of correlation is 
independent of the "reference" attenuation curve, be it 
the starburst or the SMC extcintion law. The reddening 
of the SED is therefore not due to age effects but is a 
consequence of another parameter. 

It appears that, contrary t o what has been found 
studying integrated galaxies bv lKong et alJ (12004 1 ) for in- 
stance - but in agreement with Uohnson et all (|2007l ) - 
there is no obvious trend between the perpendicular dis- 
tance and the age in the case of individual star forming 
regions. 

9.2. Necessity of a range of extinction laws 

In this study we have t ested 3 different kinds of 
extinction laws: st arburst (jCalzetti et al. | [1994 120001; 
Meurer et all 11999ft . foreground SMC (| Gordon et all 
20031) and mixed SMC. Conventionally the perpendic- 
ular distance is derived from the IRX-/3 s tarburst rela- 
tion (jKong et al.ll2004HCortese et al.ll2006l . for instance). 
There are several indications that there is a range of ex- 
tinction laws at work in quiescent star forming galaxies. 

As inferred from Figures \7\ [T2] and [T31 both a star- 
burst attenuation curve and an SMC extinction curve 
with foreground dust are required to span the locus oc- 
cupied by the data points. Neither curve individually is a 
perfect match to the data, and both (plus combinations 
of the two) are required to span a large range of the data 
locus. Likely this indicates a variation in the applicable 
dust geometries to the individual regions. 

In section [3 we have shown that assuming a range of 
models (starburst, foreground SMC assuming Ay = Ay" 8 
or Ay = 0.44 x A g y S - this ratio is prominently depen- 
dent on the dust-stars geometry - and mixed SMC), it is 
possible to reproduce simultaneously the UV color and 
the Ha emission of most star-forming galaxies. We as- 
sume that the starburst and foreground SMC are the two 
most "extreme" laws at work in quiescent star forming 
galaxies and that there is a range of dust-stars geome- 
tries in between. The two curves (and the SMC needs 
to be implemented with differential extinction between 
stellar and nebular emission) bracket most of the range 
of observed values. A variation of the extinction curve 
ha s also been detect ed for UV-luminous galaxies at z ~ 2 
bv lNoll et all (l200l . 

9.3. Influence of the star formation history and dilution 
by an older population 

Dilution of actively (Ha emitting) star forming regions 
by an older (~ 300 x 10 6 years) instantaneous burst pop- 
ulation is a necessary requirement for a few of the reddest 
(in UV colors) regions but not for all. 

We remove local background from all our star form- 
ing regions, thus removing Hubblc-timc-avcraged stel- 
lar populations. However, the typical region sampled 



by our apertures is around 300-500 pc. In this physical 
region, we can expect multiple populations to be present 
and possibly contribute to the observed colors. Indeed 
for the reddest UV color regions, the required mass ra- 
tio between an ionizing stellar population and an older, 
300 x 10 6 years old stellar population is about 1/5. This 
indicates the possibility of 5-6 generations of stars over 
the past 300 x 10 6 years, within a 500 pc region. However 
our data cannot discriminate between this possibility and 
the alternative of yet more complex star-dust geometries 
and extinctions than what considered here. 

Interestingly we note that very few regions that are de- 
tected in the UV do not have any Ha counterpart. One 
of the possible reasons is that regions that are sufficiently 
old do not emit in Ha anymore and tend to be dissolve to 
become part of the diffuse UV emission from the galaxy. 
Hence they were not selected and only acted as a contam- 
inant. The fact that we do not find the trend that can be 
found in integrated galaxies hints that it may be due in 
reality due to the contamination by older stellar popula- 
tions. Indeed, as we have seen in this paper the age or 
the latest generation of stars does not seem to play any 
significant role. We think the birthrate parameter is not 
appropriate either in the sense that it is inaccurate as it 
averages the SFR over the life of the galaxy whereas the 
typical lifespan of the UV is about 10 9 years. We there- 
fore think that the fact that few regions are not detected 
in Ha hints in favor of the scenario of the contamination. 

Finally, it can be noted that Dale et al. (2009) found 
a trend between the perpendicular distance and the Ha 
equivalent width calculated using the same method as 
the one in the paper. The most likely explanation is that 
the trend is due to the pollution of the underlying con- 
tinuum by evolved stellar populations. Indeed, their av- 
erage Ha equivalent width are averaged over the galaxy. 
If we assume that the mean age of star forming regions 
is the same in all galaxies - which is is a reasonable as- 
sumption for quiescent star forming galaxies -, the the 
differing amount of underlying stellar populations can 
create such a correlation: they simultaneously increase 
the perpendicular distance reddening the UV color and 
lower the Ha equivalent width by increasing the flux of 
the continuum. In addition, as the R-band emitting and 
the UV emitting stars have a different life expectancy, it 
induces some scatter in the relation. 

9.4. Caveats of star formation rate derivations 

From what we have shown previously, the usual star 
formation rates estimates based on the assumption of 
the starburst IRX-/3 relation are unreliable in the case 
of quiescent star forming galaxies. To estimate the star 
formation rate one approach can be to derive an empir- 
ic al law from a s ample of galaxies, like has been done 
bv lCortese et al.l ([2006D for instance, however the scatter 
remains important. 

To determine accurately the actual law, it is crucial 
to estimate the fraction of the UV color which is con- 
tributed by an older population as we have shown in sec- 
tion 16.51 Taking into accoun t the age of the stellar pop- 
ulation, ICortese et al.l (|2008f ) for instance have provided 
some recipes to correct for the UV extinction, assuming 
a Large Magellanic Cloud extinction law. Another crit- 
ical point to derive the actual star formation rate is the 
choice of an extinction law. Indeed, as we have already 
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mentioned, the effect of the extinction is much more dra- 
matic on the UV color in the case of a SMC law. This 
can be seen in Figure [12] for instance. In Figure [16] we 
plot the ratio of the estimated SFR assuming a starburst 
attenuation law to the estimated SFR assuming a SMC 
extinction law as a function of the UV color. 




0.4 0.2 0.0 -0.2 

log(F x (FUV)/F x (NUV)) 



Fig. 16. — Ratio between the SFR derived from the IRX-/3 re- 
lation for the starburst attenuation law (SFR(SB)) and for the 
foreground SMC extinction law (SFR(SMC)) as a function of the 
UV color and the extinction (in magnitudes) of the stars. The color 
indicates the value of the ratio. Assuming an extinction lower than 
1.5 magnitude, SFR(SB) can overestimate the SFR by a factor 2 
to 10.5 relatively to SFR(SMC). The white lines indicate the value 
of the ratio. The black lines indicate the track for a given age. 

We see that indeed the effect on the derived star forma- 
tion law can be quite dramatic. If we assume a starburst 
attenuation law, for Ay = 1 we can overestimate the 
SFR by a factor up ~ 7 if the actual extrinction law is 
foreground SMC like. Conversely, a similar error, but in 
the sense of underestimating the SFR, will be made if 
the foreground SMC extinction is adopted, but the ac- 
tual extinction is described by a starburst attenuation 
curve. 

As we have shown, the age seems to be only a mi- 
nor factor regarding the location on the IRX-/3 diagram. 
However, standard star formation rate estimators a re sig - 
nificantly age-dependent. Indeed, the iKennicuttl (]1998f ) 
estimators for instance are based on the assumption of 
a constant star formation rate during 100 x 10 6 years. 
However if it is sensible for a galaxy-averaged luminosity, 
it is unlikely that small individual star forming regions 
undergo such a star formation history. If we assume a 
foreground SMC, the best-fit age is 6 x 10 6 years. 

9.5. Completeness of the parameter space sampling 

Even if we have selected over 300 star forming regions, 
the location of each of the region is dependent on the 
galaxy metallicity for instance. As no such study has 
been performed before, we can only compare to samples 
of integrated galaxies. In Figure [7] we plot the IRX-/3 
diagrams along with data points from SINGS and LVL 
(Dale et al. 2009). 

We see that the IRX range is very similar to that of 
individual star forming regions in galaxies. Regarding 
the UV color, it tends to be redder. This is particularly 
striking for galaxies which have a small IRX lower than 



—0.3, but have red UV colors of 0.24 ± 0.07 compared 
to 0.37 ± 0.05 for galactic star forming regions. This can 
be explained by the influence of an older stellar popula- 
tion. Indeed, photometry of integrated galaxies is par- 
ticularly influenced by a near-ultraviolet emitted older 
stellar population which reddens the UV color. Con- 
versely, individual star forming regions are relatively less 
affected by older populations as 1) the specific star for- 
mation rate in a star forming region is higher than for an 
integrated galaxy hence lowering the relative luminosity 
of any older population 2) the way the background level is 
measured naturally subtracts the contribution of an older 
populations in the aperture (see section [?TT])- Therefore, 
rather than an undersampling of the parameter space, 
the fact that we have not detected regions which have 
simultaneously a red UV color and a low IRX is a con- 
sequence of a bias introduced by the emission of older 
stellar populations in the case of integrated galaxies. As 
we can see in Figure the envelope of SINGS and LVL 
galaxies is shifted towards redder colors compared to the 
individual star forming regions and there is a significant 
number of galaxies with extremely red UV colors: 9.6% 
of the galaxies in the combined LVL and SINGS sample 
and 1.2% of the individual star forming regions have a 
log F\ (FUV) /F x (NUV) redder than 0. 

10. SUMMARY AND CONCLUSIONS 

We have performed an analysis of over 300 star form- 
ing regions in 8 local quiescent star forming galaxies 
to garner new insights on the IRX-/3 relation by using 
Spitzer infrared, GALEX ultraviolet and ground based 
optical broad band and narrow-band Ha data. Specifi- 
cally we investigate the deviations of star forming regions 
in normal star forming galaxies from the locus identified 
by the starburst attenuation curve, in hope of obtaining 
information on similar deviations by normal star forming 
galaxies. 

The mean age of the stellar population, as traced by 
the U — B color and by the equivalent width of Ha, 
is shown to be relatively uncorrelated or only weakly 
correlated with the perpendicular distance of the data 
from the starburst attenuation curve in the IRX-beta 
diagram. Stellar population ages have been suggested 
as a potential culprit for the dev iations in the cas e of 
normal star-forming galaxies fe.g. iKong et al.ll2004l . via 
the b-parameter), but we do not find evidence for such 
trend in HII knots of galaxies. Our results are in agree- 
me nt with similar resu lts, o btained for wh o le gal axies, 
by iSeibert et all (|2005Q and Uohnson et all (|2007f ). but 
disagree with recent findings by Dale et al. (2009), also 
derived for whole galaxies. 

We find that a range of dust extinction/attenuation 
curves and dust/star geometries are actually required to 
account for the location on the IRX-beta plot of most of 
our HII knots. Specifically, we find that the starburst at- 
tenuation curve and the SMC extinction with foreground 
geometry bracket the location of most of the data points, 
not only on the IRX-beta diagram but on other plots re- 
porting colors and other properties of the regions. For 
the reddest (in UV color) regions, we need to include 
contributions to the observed colors and fluxes of under- 
lying older (~ 300 Myr) stellar populations to account 
for the observed UV and U — B colors. The mass ratio 
between the young (ionizing) stellar population and the 
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older one is about 1/5. 

Individual star forming regions and integrated galaxies 
do not populate the same locus in the IRX-/3 diagram. 
Indeed, the former have a bluer UV color and are more 
tightly correlated. The most likely explanation is that 
older stellar populations up to 1 - 2 x 10 9 years old, 
which have a redder UV color redden the observed UV 
color of star forming galaxies. To accurately determine 
the star formation rate, it is crucial to take into account 
this older population. 



Finally, a strong correlation is found with the metal- 
licity, in the sense that the more metal rich regions tend 
to display redder UV colors and larger IR/UV values, 
as expected metallicity and dust content correlate with 
each other. 



This work has been supported by NASA ADP grant 
NNX07AN90G. 
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APPENDIX 
APERTURES 



In Figure [T71 we see the apertures ovcrplotted on a section of NGC 3031. This illustrates the fact that some regions 
are selected in FUV only or 24 /im only. Some regions though are selected from both bands but the shape of the 
aperture may vary depending on the actual shape of the star-forming region in a given band. 



* 



.ii 



Fig. 17.— Section of NGC 3031 seen in FUV (blue), Ho (green) and 24 /im (red), 
from FUV (green) and 24 fim (red). 



The overplotted countours are the apertures made 
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